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MARTIAN SEISMICITY 
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Introduction. The design and ultimate success of network seismology experiments on 
Mars depends on the present level of martian seismicity. Volcanic and tectonic landforms 
observed from imaging experiments show that Mars must have been a seismically active 
planet in the past and there is no reason to discount the notion that Mars is seismically active 
today but at a lower level of activity. Here we explore models for present-day seismicity. 

Method. The frequency of occurrence, N f , of marsquakes as a function of surface wave 
magnitude, M s , is given by the relationship 

dNf/dMj. = Afexp(-BM S ) (1) 

where the B-value adopted is 2.07 [1], The maximum rate of shear strain over a time inter- 
val t can be related to the sum of the seismic moments, EMq, through [2] 

d^max/dt = EMo/pVt (2) 

where V is the seismogenic volume and \i is the shear modulus (here taken at 7xl0 10 Pa). 
The sum of the seismic moments can be expressed as 
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IMo= Mo(Ms)[dN^dMs]dMs 
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where Nj is the rate of occurrence of marsquakes (a function of Ms) and Mg(Max) is the 
magnitude of the largest marsquake, which we guess by comparison to the Earth and Moon. 
We use a relationship between seismic moment and shear wave magnitude for intraplate 
earthquakes in the terrestrial oceanic lithosphere [3] (these events are probably thermoelastic 
in origin): 


M 0 (M S ) = 4.68X10 1 1 exp(2.58M s ) (4) 

where the units are N-m. Substituting (1), (2), and (4) into (3) yields in units of (yr- 
magnitude)'! 


A f = 2.40x1 0 6 [(de max /dt)V]exp[-0.509M s (Max)] (5) 

Substitution of (5) into (1) and integrating over the interval (Mg, Ms + AM § ) yields the rate 
of occurrence of marsquakes in that interval as a function of M s (Max), strain rate, and seis- 
mogenic volume. We proceed to estimate strain rates and the attendant seismogenic 
volumes for various processes operating in the martian environment. 

Sources of Global Seismicity. Sources of martian seismicity include thermoelastic 
cooling, regional loading of the surface including Tharsis and the polar caps, changes in the 
principal moments of inertia, obliquity changes, daily and annual solar tides, and atmo- 
spheric coupling. Seismicity is estimated for the first three mechanisms. 

We performed a simple parameterized convection calculation [4], which shows that the 
mantle of Mars has cooled about 3° C and the lithosphere has cooled about 10° C (and thick- 
ened) in the last 10 8 years. In the oceanic lithosphere, earthquakes occur as deep as the 800° 
C isotherm [5]. In the martian thermal model this corresponds to a depth of about 150 km, 
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which defines the seismogenic thickness, Z L . The cooling rate in the lithosphere can be re- 
lated to strain rate through the coefficient of thermal expansion, a. If Z L « R<), the martian 
radius, then 


d^max/dt - a(dT/dt)Z L /R 0 (6) 

Using NUMax) = 6, strain rate as calculated from equation (6) (= -6x1 O' 21 s’ 1 ), and a 
spherical shell of thickness Z L for seismogenic volume yields the frequency distribution for 
marsquakes shown in the figure below. From this distribution we see that thermoelastic 
cooling of the lithosphere would be expected to yield 17 events greater than magnitude 4 in 
a decade. This result is probably a lower bound in that much higher cooling rates without a 
proportional decrease in seismogenic volume might be expected in the Tharsis region. 
Whether such events are detected in a network depends, inter alia, on the seismic attenua- 
tion and scattering properties of the martian interior [6]. . 

If the Tharsis load is supported in part by thermal buoyancy associated with a magmatic 
source region, then loss of buoyancy due to cooling leads to a strain rate: 

de/dt * -pa(dT/dt)g 0 h/E O) 


where go is planetary gravity, h is the vertical dimension of the cooling region, and IE is 
Young’s modulus. Using h = 200 km and the same cooling rate as above leads to dE/dt = - 
lxlO* 21 with E = 2x10* 1 Pa. Combined with a smaller seismogenic volume, this process ap- 
pears to be less important than thermoelastic cooling. The cooling and contraction of 
Tharsis will also change the difference in the principal moments of inertia, but the ac- 
companying strain rate is also smallJrelative to thermoelastic effects per se. 

Conclusions. Depending on the sensitivity and geometry of a seismic network and the 
attenuation and scattering properties of the interior, it appears that a reasonable number of 
martian seismic events would be detected over the period of a decade. The thermoelastic 
cooling mechanism as estimated is surely a lower bound, and a more refined estimate would 
take into account specifically the regional cooling of Tharsis and lead to a higher frequency 
of seismic events. This analysis should be repeated using the moment-magnitude rela- 
tionship for body waves, and other mechanisms, as listed above, should also be evaluated. 


References. [1] F.D. 
Stacey, Physics of the 
Earth, 113 (John Wiley, 
New York, 1977); [2] S.R. 
Bratt, E.A. Bergman, and 
S.C. Solomon, /. Geophys. 
Res., 90, 10249 (1985); 
[3] E.A. Bergman, Tec- 
tonophysics, 132, 1 

(1986); [4] R.J. Phillips 
and M.C. Malin, in Venus, 
159-214 (University of 
Arizona Press, Tucson, 
1983); [5] D.A. Wiens and 
S. Stein, J. Geophys. Res., 
90, 6455 (1983); [6] N.R. 
Goins and A.R. 
Lazarewicz, Geophys. 
Res. Lett., 6, 368 (1979). 


FREQUENCY DISTRIBUTION 
THERMOELASTIC MARSQUAKES 



14 



